Dynamics of quasisolitons in degenerate fermionic gases 
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We investigate the dynamics of the system of multiple bright and dark quasisolitons generated in 
a one-component ultracold Fermi gas via the phase imprinting technique in terms of atomic orbitals 
approach. In particular, we analyze the collision between two bright quasisolitons and find that 
quasisolitons are subject to the superposition principle. 



Solitons are already well established phenomena in the 
area of dilute atomic quantum gases. They have been 
observed in several experiments with trapped bosonic 
atoms as well as described in many theoretical papers. 
The latter usually involves the mean-field approach and 
the resulting Gross-Pitaevskii equation for the conden- 
sate wave function. The intrinsic nonlinearity of this 
equation is originating from the interparticle interactions. 
As the effective interaction between atoms can be both 
repulsive or attractive, two types of solitonic excitations 
[il are possible in Bose-Einstein condensates. 

Dark solitons were realized experimentally first 
whereas there have been problems with generating bright 
solitons implied by the occurrence of the collapse in the 
attractive condensates. Finally, two ways turned out to 
be successful. In experiments of Refs. |^ |^ a large 
condensate with positive scattering length is prepared 
and next the interatomic interactions are changed to at- 
tractive ones (with small negative value of the scatter- 
ing length) by using Feshbach resonances. In this case 
the single bright soliton or a train of solitons were ob- 
served. Another way of creating bright solitons was 
demonstrated in the experiment of Ref. where the 
condensate still remains repulsive but appropriate en- 
gineering of atoms in a periodic potential allows for a 
change of the sign of the effective mass and results in 
the generation of bright gap solitons. Recently, another 
concept of bright matter waves in a repulsive condensate 
has been developed that involves degenerate Fermi gas 
as a "stabilizing medium" l8|. 

In a series of papers [lol | we have shown that sim- 
ilar structures can also be created in a one-component 
ultracold fermionic gases even though to a very good 
approximation spin-polarized fermions at zero temper- 
ature do not interact. These structures exist both in the 
systems with uniform and nonuniform densities under a 
quasi-one-dimensional confinement. They can be gener- 
ated by using the technique of phase imprinting and its 
appearance is related to the Fermi statistics. There ex- 
ists the regime of parameters, defined by the condition 
that the width of the imprinted phase is bigger than the 
Fermi length, where only states with momentum close to 
the Fermi momentum are excited. Hence, the life time 
of generated structures can be long enough (in compar- 
ison with the Fermi time) allowing for observation and 
making them similar to solitons. 



The method of phase imprinting has been already 
employed in the experiments where dark solitons have 
been generated 0, Q ■ The strong and short enough off- 
resonance laser pulse is passing first through the appro- 
priately tailored absorption plate and next through the 
sample of atoms. Under such conditions the motion of 
atoms is frozen and the only result of the interaction 
of light with the atoms is imprinting the phase on the 
atomic wave functions. After the light is gone, the den- 
sity and the phase change and desired structures (like 
solitons or vortices) appear in the system, according to 
the pattern written in the absorption plate. 

In this paper we further investigate the motion of 
fermionic quasisolitons by analytical calculations based 
on the propagator technique. These calculations allow 
us for a detailed description of the collision between 
fermionic quasisolitons. We come about the conclusion 
that quasisolitons fulfill the superposition principle, i.e. 
the sum of two quasisolitons is a two-quasisoliton solu- 
tion. 

So, we consider a system of N fermionic atoms confined 
in a one-dimensional box with periodic boundary condi- 
tions at zero temperature. We assume that the many- 
body wave function of such a system is well described by 
the Slater determinant 
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The gas is spin-polarized and therefore the s-wave scat- 
tering, dominant at low temperatures, is absent and effec- 
tively the gas becomes the system of noninteracting parti- 
cles. The many-body Schrodinger equation is then equiv- 
alent to the set of one-particle equations, each of them 
describing the evolution of a particular single-particle or- 
bital ipi{x), (fNix) under the same Hamiltonian. Thus 
these orbitals remain orthogonal during the evolution and 
the one-particle density matrix is at any time given by 
the following formula 



pi{x',x",t) 
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N 
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and its diagonal part is the particle density 



1 ^ 



(3) 
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FIG. 1: Phase pattern leading to the generation of a single 
bright-dark quasisolitons pair. 

At zero temperature particles occupy the lowest pos- 
sible levels. The corresponding eigenstates are just the 
plane waves (pn{x) = exp {i pnX / h) / y/L with momenta 
Pn quantized according to Pn — 2T:hn/L, where n = 
0, ±1,±2,.... The evolution of the single-particle wave 
function can be calculated based on the propagator tech- 
nique 



(Pn{x,t) ^ K{x,x' ,t) (pn{x' ,0) dx' , 



(4) 



where K{x,x',t) is the propagator function for a one- 
dimensional box with periodic boundary conditions 



K{x,x',t) = 



exp 



(5) 



It is convenient now to increase simultaneously the num- 
ber of particles and the length of the box in such a way 
that their ratio (i.e. the Fermi momentum) remains con- 
stant. Increasing both quantities to infinity defines the, 
so called, thermodynamic limit. In this limit the prop- 
agator takes the same form as for a free particle and is 
given by 



K{x,x',t) 




exp 



2ht 



(6) 



To simplify further calculations we assume the follow- 
ing form of the imprinted phase (see Fig. ^| 



, X < —a 

<^W = <! T(f + 1) - -a<x<a 
w , X > a 



(7) 



where the parameters a and w define the width and 
the jump of the imprinted phase, respectively. The 



Fermi momentum in one-dimensional space, defined as 
hN/ {2L), determines the characteristic length Xp via the 
relation Xp = h/pp = 2L/N. From now on we choose 
the Fermi length as a unit of length. This quantity does 
not change while realizing the thermodynamic limit pro- 
cedure. Simultaneously, the expression h/ep gives a unit 
of time. Hence, according to the formula © the single 
particle wave function evolves as 



iPnix^t) 



1 



e 2t 



dx' , (8) 
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where the wave number kn = Pnl^ was introduced 
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FIG. 2: Evolution of a single bright-dark quasisolitons pair. 
Successive frames show the fermionic density (normalized to 
N/2) of 501 atoms after writing a single phase step of height 
TT and the width a = 4Af. The corresponding times are 3, 5, 
7, and 8 in units of h/ep- 

Calculating the integral (|HJl one obtains an explicit for- 
mula for the evolution of the single particle wave function 
(in units of 1/ VXf) 




iPn{.X,t) = 



'i2L 



{(l + z)(l + e™)/2 



-^SiuLJ]}, (9) 

where C{z) and S{z) are the Fresnel integrals and 
their arguments are the following functions of position 
and time 



{a — X + kn — 

TT 

(—a — X + kn 



u„ = u„ 



W r 
27ra 

W r- 

U-a + — vt 
27ra 



(10) 
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The fermionic density (the diagonal part of one- 
particle density matrix) is just the sum of single-particle 
densities obtained from single-particle orbitals given by 
©• Fig. [21shows such a density for a system of 501 atoms 
at various times after writing a phase step of height tt and 
the width 4A_F. As already reported in Ref. 0,^3 the 
bright and the dark quasisolitons are generated. They 
propagate with different velocities and in opposite direc- 
tion. The speed of bright quasisoliton is slightly higher 
than the speed of sound (equal to pp/iTi) whereas oppo- 
site is true for the dark quasisoliton. We confirm all the 
details of dynamics of such a system described earlier in 
Ref. 0,01, especially its dependence on the ratio a/Xp 
of the width of the imprinted phase step and the Fermi 
length. 
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FIG. 3: Phase pattern used for generating two bright-dark 
quasisolitons pairs. In this case bright quasisolitons move 
towards each other and collide. 

The purpose of the present paper is, however, to in- 
vestigate the collision between bright fermionic quasisoli- 
tons. To this end, we imprint the phase pattern as in Fig. 
|31 on the fermionic gas. It results in a creation of two 
pairs of bright-dark quasisolitons with bright quasisoli- 
tons moving towards each other. It turns out that the 
time evolution of individual orbitals can be expressed in 
terms of the evolution corresponding to the phase pattern 
already considered in Fig. ^and given by the formula 
First, let us notice that the phase shown in Fig. Ocan 
be split in the following way (f)^"^ ^ (f)^ + (ff — where 
and are the left and the right cuts of the phase shown 
in Fig. 13 respectively. Next, using the formula ^ one 
gets 



(11) 



where 



A = 



1 



Vi2Li 
l + i 



gi^(3:-2;')^g»fc„a;'gito^^/ 



(12) 



and superscripts 7' and V indicate that the wave func- 
tions are propagated after the imprinting the single step 



patterns (fi and (jf , respectively. The superscript 7r' 
means the evolution according to the phase pattern plot- 
ted in Fig. El 
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FIG. 4: The relative (with respect to 1/ L) oscillations of the 
expression B + 2/L as a function of position for a particular 
(n = 200) single-particle orbital. The length of the box L — 
250Af and the time is 12 in units of h/sp- Other parameters: 
c = 15, di = 4, d2 = 5, in units of Af, and w = n (the case of 
the last frame of Fig. |SJ|. 

Both orbitals ipl^{x,t) and (p^{x,t) can be calculated 
based on the following rules that connect the time evo- 
lution due to the phase pattern given in Fig. ^with the 
evolution resulting from patterns obtained by the shift 
{ipn\x,t)), the reflection {ipn'^^\x,t)), and a combina- 
tion of the reflection and the shift {(fin^'' (^7^)) of that one 
shown in Fig. ^ 

^i'\x,t) = e^'="V„(a;-s,t) 
^(:^/)(a;,i) = ^_„(-x,t) 
^l:'Hx,t) = e-''="V-»(-a;-s,t)- (13) 

One has 

(14) 

The relation Ullfl can be rewritten now in terms of the 
single-particle densities 

Pnix, t) = p\,{x, t) + pI{x, t) + ^+B 

B = 2i?e(^U^,iX*(^,i) " ^n{x,t)A* - ^;(x,t)A*) 

(15) 



and it turns out that approximately 



Pnix.t) = Pn{x,t) + Pl{x,t) 



(16) 
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FIG. 5: Two bright quasisolitons collision. Successive frames 
show the fermionic density (normalized to N/2) of 501 atoms 
after writing two phase step pattern (as in Fig. |^ of height n 
and the widths equal to 4Af (the left one) and 5Xf (the right 
one). The corresponding times are 5, 7, 8.3, 9, 10, and 12 in 
units of h/ep- Other parameters: c = 15 in units of \f- 



In fact, what is left, i.e. the expression B + 2/L, oscil- 
lates as a function of position around zero with a relative 
(in comparison with amplitude usually of the order 
of 10^"* right after the quasisolitons are created. This 
amplitude increases in time but remains at the low level 
(10~^) even after the collision of qusisolitons. An exam- 
ple of that is given in Fig. ^ where we consider the box 
of length L — 250Aj7 (it corresponds to the last frame of 
Fig. [SJ. Only when the imprinted phase step gets very 
sharp the relative amplitude becomes as large as 10^^. 
The total density is the sum of single-particle densities 
and equals 

p"-(x,t) = p\x,t)+p^{x,t)^^. (17) 

Here, the density is normalized to the number of particles 
divided by 2. Therefore, in the thermodynamic limit, 
i.e. when the size of the box and the number of atoms 
increase to infinity while keeping their ratio constant, one 
gets the following rule 



p"-(x,i) = p\x,t)+p-{x,t)-l 



(18) 



any perturbation of the initial density. Secondly, it hap- 
pens that going from H16|l to 1)17(1 the " noise" present for 
each orbital (see Fig. 0J, even if it is large, adds destruc- 
tively causing that the relation H18I) is satisfied almost 
perfectly. It has to be noticed that the same property 
of additivity ((18(1 is shared by the solutions of the wave 
equation. Here, however, single-particle densities satisfy 
the nonlinear hydrodynamic equations according to the 
Madelung representation of the Schrodinger equation ^2 



dVn 

dt 



d_ 

dx 



= 

ft2 



1 ^2 



2 2m? J]j^ dx"^ 



p;: =0, (19) 



First, this formula generalizes the results already ob- 
tained in Rcf. ^3 (and expressed by solution (24)) to 



where p„ and w„ are the density and the velocity fields, 
respectively, associated with the orbital '^'■^{x.t). 

In Fig. 13 we illustrate the collision of two bright 
fermionic quasisolitons. Successive frames show the to- 
tal density calculated based on the formulas l(ll() . 1(14(1 . 
and ©. No interference, as a consequence of ((18() . is 
observed when both peaks meet at the center of the box 
what distinguishes this case from the case of collision of 
usual wave packets. It is not surprising, since building the 
quasisoliton we add probabilities (see ^) rather than the 
amplitudes. After the collision both quasisolitons reap- 
pear, however no phase shift, usually expected when two 
solitons collide, is observed. This is another consequence 
of formula ((18(1 . So, what we have is a solution (in a sense 
of (Q) of a set of nonlinear hydrodynamic equations l(19() 
that is a sum of quasisolitons not only well before and 
after the collision but all the time. 

In conclusion, we have investigated the dynamics of 
fermionic quasisolitons generated in ultracold fermionic 
gas by using the technique of phase imprinting. We ob- 
tained an analytic formula that describes the propagation 
of a bright-dark quasisoliton pair. Based on this formula 
we were able to study the collision of two bright fermionic 
qusisolitons and found that the quasisolitons are subject 
to the superposition principle - the sum of two quasisoli- 
tons is a two-quasisoliton solution. 
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